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Abstract 

Using the nonequilibrium Green function, we show that microwave irradiation 
can reverse the supercurrent flowing through a superconductor / quantum- 
dot / superconductor structure. In contrast with the conventional sideband 
effect in normal-metal / quantum-dot / normal-metal junctions, the photon- 
assisted structures appear near Eq = ^hw (n = ±1, ±2 • • •), where Eq is the 
resonant energy level of the quantum dot and u is the frequency of microwave 
field. Each photon-assisted structure is composed of a negative and a positive 
peak, with an abrupt jump from the negative peak to the positive peak around 
Eq = r;hu}. The microwave- induced 7r-j unction transition is interpreted in the 
picture of photon-assisted Andreev bound states, which are formed due to 
multiple photon- assisted Andreev reflection between the two superconductors. 
Moreover, the main resonance located at Eq = can also be reversed with 



proper microwave strength and frequency. 
PACS numbers: 74.50.+r, 73.63.Kv, 72.23.Ad. 

When two superconductors are weakly linked, dc current can flow even without bias 
voltage. The driven force of the supercurrent is the phase gradient in the macroscopic wave 
function of Cooper pair condensate. The supercurrent and the phase difference across the 
junction has the relation I = J c sin0, with J c > being the critical supercurrent. If the 
link area of the junction is controlled by some external conditions, J c may be enhanced , 
suppressed, or even reversed. The reverse sign of I c is referred as to tc -junction transition, 
since the minus sign can be absorbed into <fi as an internal phase shift of tc. The ground 
state of 7r-junction is <fi = tc rather than = as in usual 0-junction. This is easily seen 
from the relation I — and therefore F = —^Ic cos <fi + Fo, where F is the free energy of 
the junction. The studies of 7r-junction have not only academic interests, but also potential 
applications, e.g., realization of qubit in a superconducting loop with 0- and tc- junctions 

In addition to the intrinsic tc phase shift in the order parameter of high T c superconduc- 
tors due to the d-wave symmetry , several mechanisms were also proposed and investigated 
for carrying out 7r-junction with the conventional BCS superconductors. (1) Coupling two 
superconductors through a ferromagnetic layer. The idea can be traced back to the original 
works of Fulde and Ferrel || and Larkin and Ovchinnikov M , who independently predicted 
that superconducting order parameter can be modulated by an exchange field, and contains 
nodes where the internal phase shift is tc. With the growing interests in the hybrid system 
of superconducting / ferromagnetic (S/F) materials, 7r-junction behavior in SFS sandwiches 
and SF superlattices were intensively studied Very recently, experiments on Joseph- 

son junction inserted with a weakly ferromagnetic interlayer did observe the 7r-junction 
transition ||. (2) Coupling two superconductors by an Anderson impurity or a quantum 
dot: The early papers of Glazman and Matveev |T(J and Spivak and Kivelson showed 
that when the impurity is singly occupied due to Coulomb repulsion, the sign of Josephson 
current is opposite to that without the repulsion.. In the last decade, a number of theoret- 
ical papers were devoted to this issue fT2"| - |I5|| . Experimentally, the technique of fabricating 
Josephson junctions containing nanoparticles was available [^6j, yet no relevant results on 
7r-junction transition was reported. (3) Introducing a nonequilibrium distribution in the cen- 
tral mesoscopic region: In a mesoscopic superconductor / normal-metal / superconductor 
(SNS) structure, the quasiparticle distribution can be driven far from the equilibrium by a 
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control voltage across the N region. When the control voltage exceeds a certain value, the 
nonequilibrium distribution has so much weight on the negative part of the current carrying 
density of states that the supercurrent reverses its sign. After a few years of the prediction 
fi7H20| , the reverse of the supercurrent was successfully observed in a controllable Josephson 
junction pif . 

In this paper, we propose a new mechanism for the realization of 7r-junction, schemat- 
ically shown in the inset of Fig.l. Consider a Josephson junction consist of two BCS su- 
perconductors (S) coupled through a quantum dot (QD). Applying a microwave (MW) on 
the QD, the resonant energy level of QD will shift adiabatically with the time dependent 
external field. We shall show below that 7r-junction transition occurs in the photon-assisted 
Josephson current, and the main resonance can also be reversed by a proper choice of MW 
strength and frequency. The MW-induced 7r-j unction transition is related to the formation 
of photon-assisted Andreev bound states (PAABS), which is a generalization of the usual 
Andreev bound states (ABS). 

We model the S-QD-S system by the following Hamiltonian, 



H 



Hl + H D + H R + H T 

A /3fcT u /3A : 



with (3 



(1) 

L I R is the 



where Hp 

standard BCS Hamiltonian for the left / right superconducting electrode; Ho = (E + 
W cos uot) c\c a is for the QD under MW irradiation, in which the intradot interaction 
is ignored for simplicity f52]] , and H T = p ka (v 'pap ka c a + H.c.) is the tunnel Hamiltonian , 
connecting the three parts together. 

The time dependent current flowing through S-QD-S can be formulated using the 
nonequilibrium Green function 



2 Re jrr a z j dt x [G r (t, tjE^fat) + G K {t, ti)£^i, t 
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in which the full Green function of QD and the self energy induced by the coupling with /3 
electrode are defined as 
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and cx z in Eq.(2) is the third Pauli matrix, with the diagonal element +1 for the electron 
current (spinf) and -1 for the hole current (spin|). G r and G < obey the corresponding 
Dyson equation and Keldysh equation: 

(6) 
(7) 



G r (t 1 ,t 2 ) = g r (t 1 ,t 2 ) + j j dtdtg r (t 1 ,t)J: r (t,t)G r (t,t 2 ) , 

G < (t 1 ,t 2 ) = J J ' dtdt'crfat^fat'Wit' \t 2 ) . 



The remaining task is to solve these equations properly. It should be pointed out that finite 
perturbation expansion of Eq.(6) is inadequate in the problem, because the formation of 
PAABS involves up to infinite order of tunneling processes. 

Although G r (ti,t 2 ) is no longer the function of t± — t 2: it still holds that G r (ti + ^, t 2 + 
^) = G r {ti,t 2 ) due to the periodical time dependence of the MW. Hence G r (ti,t 2 ) can be 
Fourier expanded as 



G r {t u t 2 ) = $> u <^ J ^e-^-^GKe) . 



(8) 



Define the Fourier transformed Green function 



G„„(e) = G r m _ n (e - nu) , 

which has the property that if C(t 1 ,t 2 ) = J dtA(t 1: t)B(t : t 2 ) then C mn (e) 
J2 k A m fc(e)B fen (e). The Fourier transformed g r and £ r,< can be obtained as 
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(10) 
(11) 



in which J„(a:) is the nth Bessel function with the argument a = S < (e) = /(e)[S a (e) — 
S r (e)] with /(e) being the Fermi function, and S r (e) = S£(e) + E£(e) with 
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In the self energy, r p is the coupling strength between QD and (3 electrode, iq is the dephas- 
ing rate in the S electrodes which determines the broadening of PAABS. In the numerical 
calculation, the limit 77 — > + is achieved by extrapolation. 

The Fourier transformed Dyson equation can be expanded as 



G r (e) = g r (e) + g r (e)Sr(e)g r (e) + g'"(e)Sr(e)g r (e)Sr(e)g'-(e) + 



(13) 



To re-sum up the series, we adopt the resonant approximation: [24 
1 1 
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(e - l[u + E + i0+) (e - 1' 2 oj + E + i0+) ' ^ (e - l' x u + E + i0+) 2 ' 

which is justified for the weak coupling case of Tp <C u. Under this approximation, the 
Dyson equation can be exactly solved as 
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The time dependent current is Fourier transformed as 



(16) 
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With the relations G<(e) = G r (e)£<(e)G a (e), G a (e) = [G r (e)] f , and G r (e) solved in 
Eq.(16), the current can be evaluated. In this work, we concentrate on the time aver- 
aged current / = if = -if. Notice that / = (v R lf - T L lf^j / (T L + T R ), and de- 
fine £ = (r R S L — r^S/j) / (T L + T R ). Again, we adopt the resonant approximation in 
G < (e) = G r (e)S < (e)G a (e), and the current formula is simplified as 

1= \\ Ji(«)^(«)2Re Tr a z [G\,{e)Y<{e) + G< (e)£ a (e)] . (23) 

1 1' 

where G r u ,(e) = ^(5^, (e)) 1 — E^, (e) and G<,(e) = G^, (e)S^, (e)G^, (e). One can see in the 
formula that I is contributed by various photon-assisted Andreev reflections, in which an 
electron absorbs / photons and a hole emits I photons. Eq.(23) is the central result of this 
paper. 

Below we shall numerically study the influence of the MW field on the Josephson current 
flowing through S-QD-S. In the calculation, we set e = h = 1, take k R T = 0, and choose 
A = 1 as the energy scale. In a qualitative discussion, we regard I c = I((f> — |) as the critical 
supercurrent, and use the sign of I c as the criterion of the 0-junction vs the 7r-junction. Fig.l 
shows the curves of the critical supercurrent I c vs the resonant energy level Eo, for different 
MW strengths. There are two remarkable features in the plot: (1) With the increase of 
MW strength, the main resonance located at E = is gradually suppressed, while several 
photon-assisted tunneling structures grow up around E = \nt\jjj {n = ±1, ±2 • • •). (2) Each 
of the structures contains a negative and a positive peak, and I c reverses its sign abruptly 
near Eq = ^nhw. Meanwhile, the magnitudes of these structures exhibit a non-monotonous 
dependence on the MW strength. Feature (1) is quite different from the sideband effect 
in N-QD-N, where photon-assisted peaks appear near E = nhw. In N-QD-S junctions, 
however, photon-assisted peaks also appear near E = ^nhu, which can be explained in 
the picture of photon-assisted Andreev reflection (PAAR) ||24j| . Basically, feature (1) can 
be understood in the similar way as in N-QD-S, but should be noticed that in S-QD-S the 
photon-assisted electrons and holes undergo up to infinite order of PAAR processes, so that 
PAABS are formed in the QD (schematically shown in the inset of Fig.l). Feature (2) is 
dramatically distinguished from those in either N-QD-N or N-QD-S, indicating that MW 
irradiation can lead to 7r-junction transition in the photon-assisted current through S-QD-S. 

To understand the feature (2) requires the knowledge of PAABS . For this purpose, 
the time averaged current is rewritten in the form I = |J ^/(e)j(e). The analysis of 
the current carrying spectrum j(e) provides the information of the supercurrent carried by 
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each of the PAABS. Fig. 2 illustrates the spectrum of j( e ) f° r MW strength a = 0.6, with 
E = 0, 0.05, 0.10, 0.15, 0.20. For comparison, the spectrum of E = and E = 0.2 
without MW irradiation are also shown on top of Fig.2. Without MW irradiation, the 
spectrum of S-QD-S has two ABS within the superconducting gap, located near e = ±E , 
carrying supercurrent with opposite signs. The spectrum also has continuous parts outside 
the superconducting gap, negative for u < —A and positive for uo > A. Turning on the MW, 
the original two ABS split into two sets of PAABS, due to various photon-assisted processes. 
The PAABS marked with A n locate near e = E + nu>, carrying supercurrent with positive 
signs; while those marked with B n locate near e = — E + nu, carrying supercurrent with 
negative signs. On the contrast, MW has negligible effect on the continuous spectrum, 
because the contribution from these parts depends strongly on the density of states in the 
S electrodes and hence photon-assisted processes is largely suppressed. Notice that at zero 
temperature I is contributed by the spectrum in the range of e < 0. For E = (Fig. 2a), 
the states A n and B n are energy degenerate, carrying opposite supercurrent weighted by a 
distribution factor similar to the sideband in N-QD-N. The contributions from A_i 

and B_ 1; A_ 2 and B_ 2 , etc. cancel with each other exactly, and the positive contribution 
from A overwhelms the negative contribution from C~, resulting in a positive current. With 
the change of the resonant energy level E , A n and B n move toward — A and +A, respectively. 
For E = Tjfkd (Fig.2d), the states A n and B n _i are energy degenerate, and therefore strongly 
hybridized with each other. So the weight of A n and B n _i are re-distributed, such that they 
carry the same amount of supercurrent with opposite signs. The contributions from A 
and B_i, A_i and B_ 2 , etc. cancel with each other exactly, leaving the negative continuous 
spectrum C~ contribute to the supercurrent. This is the origin of the 7r-junction transition in 
the photon-assisted supercurrent. Further investigation shows that the evolution of PAABS 
from (a) to (d) is more complex than stated above. In fact, the hybridization of A n and B n _! 
already occurs before the energy degeneracy (see Fig. 2d and its inset). The hybridizing of 
PAABS obeys the following rule: when two states Aj and Bj approaches each other, they 
"interact" with each other only if the energy of Aj is higher than that of Bj. With this rule, 
the abrupt sign change in the current near E = ^fku is readily understood. 

In the above, we have shown that MW irradiation can reverse the sign of photon-assisted 
Josephson current. We are still curious whether the main resonance located at E = can 
also be reversed. In the case of E = 0, the current formula is further simplified since I = t 
is required by the resonant tunneling. Fig.3 shows the curves of I = 7(0 = §,£b — 0) vs 



7 



the MW strength a, for four groups of MW frequencies. These numerical results reveal some 
of the properties of Iq. (1) Iq has non-monotonous dependence on a, which is natural due 
to the oscillatory nature of Bessel function. (2) For some frequencies, e.g., to = 0.3 in Fig.l, 
Iq is always positive; while for others, e.g., uj = 0.6, Jo reverses its sign in a certain range 
of a. In other words, MW with proper strength and frequency may also induce 7r-j unction 
transition on the main resonance. (3) The dependence of Iq on a changes abruptly around 
the frequency of uj — A/n (n — 1,2,3 • • •). Roughly, the reason is that when the frequency 
crosses A/n, A± n and B± n jump out of the superconducting gap and feel the singularity in 
the density of states of S electrodes. 

In conclusion, we have shown theoretically that MW irradiation has a non-trivial effect 
on the Josephson current flowing through S-QD-S structure.. The photon-assisted tunneling 
structures appear near E = ^ftw, and each contains a negative and a positive peak, indi- 
cating 7r-junction transition induced by MW. The main resonance located at E = may 
also reversed with proper MW strength and frequency. The understanding of these results 
requires the knowledge of PAABS, which are formed due to multiple PAAR. We hope that 
this work can stimulate experimental interests in S-QD-S structure, where QD may be a 
nanoparticle or a gate defined geometry in 2DEG. 

This project was supported by NSFC under Grant No. 10074001. T. H. Lin would also 
like to thank the support from the State Key Laboratory for Mesoscopic Physics in Peking 
University. 
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FIGURE CAPTIONS 

Fig. 1 The curves of the critical supercurrent I c vs the resonant energy level E , with the 
increase of MW strengths. Parameters are: u = 0.3, T = 0.01. The inset in a = 
schematically shows the setup of an S-QD-S structure under the MW irradiation. The 
inset in a = 0.3 demonstrates the formation of PAABS around E = —^huj. 

Fig. 2 The analysis of current carrying spectrum j(e) for the point (a)-(e) marked in the 
curves of a = 0.6 in Fig.l. For comparison, the spectrum for E = (solid) and 
E = 0.2 (dotted) without MW irradiation are also shown on the top. 

Fig. 3 The curves of the main resonance Iq vs the MW strength a, for MW frequencies uo 
chosen near A/3, A/2, A, and uj > A. 
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